Measurement of production cross section by 8 GeV protons using high energy 
neutrino interactions in the SciBooNE detector 
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The SciBooNE Collaboration reports K^ production cross section and rate measurements using 
high energy daughter muon neutrino scattering data off the SciBar polystyrene (CgHg) target in the 
SciBooNE detector. The K^ mesons are produced by 8 GeV protons striking a beryllium target in 
Fermilab Booster Neutrino Beam line (BNB). Using observed neutrino and antineutrino events in 
SciBooNE, we measure 

-r-^ = (5.34 ± 0.76) mbl{GeVlc x sr) 
dpail 

for p + Be — ^ K^ + V at mean K'^ energy of 3.9 GeV and angle (with respect to the proton 
beam direction) of 3.7 degrees, corresponding to the selected if"*" sample. Compared to Monte 
Carlo predictions using previous higher energy K^ production measurements, this measurement, 
which uses the NUANCE neutrino interaction generator, is consistent with a normalization factor of 
0.85±0.12. This agreement is evidence that the extrapolation of the higher energy K^ measurements 
to an 8 GeV beam energy using Feynman scaling is valid. This measurement reduces the error on 
the A"*" production cross section from 40% to 14%. 
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I. INTRODUCTION 

Inclusive kaon production by low-energy protons (1 to 
15 GeV) is of interest both theoretically and experimen- 
tally. In this low-energy region, kaon production is dom- 
inated by exclusive processes. For example, the lowest 
threshold K'^ production process is p + p — > + A+p, 
which for a fixed target setup has an incoming beam 
energy threshold of 2.52 GcV. Since exclusive channel 
threshold effects are jmportant, theoretical models such 
as Feynman scaling jjy] may be better in describing low- 
energy production cross sections. Measurements of kaon 
production in this region are not extensive and do not 
cover wide kinematic regions. In addition, systematic 
data on the energy and target nuclei dependence is not 
available. Thus, new measurements of kaon production 
are needed in this region. Experimentally, kaon produc- 
tion is also relevant for neutrino experiments since im- 
portant components of the incident neutrino flux come 
from kaon decays. 

A primary motivation of this work is to verify the sim- 
ulation of neutrinos from the Fermilab Booster Neutrino 
Beam (BNB) line with actual data. The BNB line pro- 
vides neutrinos for the MiniBooNE and SciBooNE [1] 
experiments, as well as possible future experiments, in- 
cluding MicroBooNE Q. In this beam line, protons with 
8 GeV kinetic energy are directed onto a 1.8- interaction 
length beryllium target. The average energy of 7r+ {K'^) 
that decay to neutrinos in the MiniBooNE detector ac- 
ceptance is 1.89 (2.66) GeV. Therefore 37.6% (92.1%) 
decay before the end of the 50 m long decay region. 
The relevant decay modes for MiniBooNE/SciBooNE are 
7r+ — >■ fi^i^fj,, — i> fJ,~^tyfj., which produce 92.9% of the 
neutrino beam, tt~ — > ijl~v^, which produces 6.5% of 
the neutrino beam, and K'^ — 7r'^e"*"i^e, /^^ — ^ g^v^v^, 
— >■ Tr~e~^Ve, and -> n^e~De, which produce 0.6% 
of the neutrino beam. 

While the neutrino flux is predominantly due to 
decay, decay is the dominant source above 2 GeV. 
The neutrinos from kaons provide a unique source of high 
energy events for experiments on the BNB line studying 
neutrino cross sections, and can represent a source of 
background for experiments exploring neutrino oscilla- 
tions and beyond-the-standard-model effects. Therefore, 
it is important for the BNB line experiments to under- 
stand the rate of production. 

An accurate understanding of production will re- 
duce systematics associated with the measured Ve back- 
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ground in MiniBooNE, a major contributor to the un- 
certainty in the previously published Ve oscillation ap- 
pearance result [5|. The measurement of produc- 
tion in this energy region combined with K"^ production 
at higher energies is a good test of production models 
such as the Feynman-Scaling [6] and Modified Sanford- 
Wang [rj parameterization used to describe secondary 
meson production at low primary proton beam energy in 
the BNB. 

This work describes a measurement of K'^ production 
by measuring the rate of high-energy events from kaon 
decay. The data sample used for the measurement comes 
from the interaction of z/^ and which undergo charged 
current (CC) neutrino interactions in the fiducial volume 
of the SciBooNE detector, generating high energy /i~ 
and fi^ (along with a host of other particles) that pene- 
trate the entire SciBooNE detector, providing essentially 
a minimum muon momentum requirement of 1.0 GeV/c. 
The neutrinos from decay can be isolated using the 
angular distribution of the outgoing muons and the mul- 
tiplicity of charged particles produced in the interaction. 
The number of decay neutrinos is then compared to 
prediction to make a determination of a normalization 
factor for production in the BNB and correspondingly a 
production cross section. The paucity of high-energy 
events in the SciBooNE experiment prevents a measure- 
ment of the kinematic distribution of production but 
does allow a normalization determination with improved 
uncertainty. 



II. SCIBOONE EXPERIMENT 
A. Neutrino Beam 

The SciBooNE experiment detected neutrinos pro- 
duced by the Fermilab BNB. The same BNB beam is also 
serving the MiniBooNE experiment. The BNB uses pro- 
tons accelerated to 8 GeV kinetic energy by the Fermilab 
Booster synchrotron. Beam properties are monitored on 
a spill-by-spill basis, and at various locations along the 
BNB line. Transverse and directional alignment of the 
beam, beam width and angular divergence, beam inten- 
sity and losses along the BNB, are measured and used 
in the data quality selection. Protons strike a 71.1 cm 
long beryllium target, producing a secondary beam of 
hadrons, mainly pions with a small fraction of kaons. A 
cylindrical horn electromagnet made of aluminum sur- 
rounds the beryllium target to sign-select and focus the 
secondary beam. For neutrino running mode, the horn 
polarity was set to focus particles with positive electric 
charge and for antineutrino running mode, the horn po- 
larity was set to focus particles with negative electric 
charge. The neutrino beam is mostly produced in the 
50 m long decay region. The analysis described in this 
paper will use data from both neutrino and antineutrino 
running modes. 

The Monte Carlo (MC) simulation of the neutrino 
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beam was modeled by the MiniBooNE collaboration. 
The MiniBooNE collaboration uses a GEANT4-based 
Monte Carlo simulation that can be roughly divided into 
five consecutive simulation steps. The first simulation 
step is the definition of the beam-line geometry includ- 
ing the shape, location, and composition of components. 
The second simulation step is the generation of primary 
protons according to the measured beam optics proper- 
ties. The third simulation step is the simulation of par- 
ticles produced by the initial p-Be interaction, including 
the elastic and quasi-elastic scattering of the protons in 
the target. Custom tables for the production of proton, 
neutron, tt^, K^, are used based on the phenomenol- 
ogy of particle production and data of the production of 
these particles at higher energies. The fourth simula- 
tion step is the propagation of the produced particles us- 
ing the GEANT4 framework taking into account energy 
loss, electromagnetic and hadronic processes, and trajec- 
tory deflection by the magnetic field created by the horn. 
The fifth simulation step is the decay of propagated par- 
ticles into neutrinos using the current branching fraction 
measurements 2]. 

Particle production is simulated using the methods de- 
scribed in Ref. [2] . The production of is simulated us- 
ing a Feynman scaling formalism based on p-Bepro- 
duction data at different primary proton energies [J, 0] . 
The predicted double differential cross section using the 
Feynman parametrization reported in 0, is 



— — = (6.3 ± 2.5) m6/(GeF/cx sr), (1) 
dpdil 

at the mean energy of 3.9 GeV and mean angle 
of 3.7 degrees, which are the mean energy and angle for 
kaons which produce neutrinos in SciBooNE. 

For 7r+ and vr^ production, the Sanford-Wang (SW) 
parametrization to the HARP p-Be data [g| at 8.89 
GeV/c is used to determine the central value with as- 
sociated errors determined from spline fits. SW pro- 
duction is also used for production and errors. The 
long life time of the X£ combined with the fact that 
they are not focused by the magnetic horn leads to the 
expectation that the contribution of decay neutrinos for 
this source is small relative to the . For produc- 
tion, the scarcity of production measurements in the rel- 
evant kinematic regions motivated the use of the MARS 
hadronic interaction package Q to determine the abso- 
lute double differential cross-sections. 

The neutrinos produced from the simulation are ex- 
trapolated along straight lines toward the SciBooNE de- 
tector. All neutrinos whose ray traces cross any part of 
the detector volume and surrounding rock are considered 
in the SciBooNE neutrino flux prediction and the kine- 
matics of the neutrino and their parents are stored. 

In the neutrino mode running (positive horn polarity), 
a total neutrino flux of 2.2 x 10~^cm~^/POT is expected 
at the SciBooNE detector location, with a mean neutrino 



energy of 0.7 GeV. The flux is dominated by muon neu- 
trinos (92.92% of total), with small contributions from 
muon antineutrinos (6.48%), electron neutrinos (0.54%) 
and electron antineutrinos (0.05%). 

In the antineutrino mode running (negative horn po- 
larity), a total neutrino flux of 1.3 x 10~^cm~2/pOT 
is expected at the SciBooNE detector location, with a 
mean neutrino energy of 0.6 GeV. The flux is dominated 
by muon antineutrinos (83.85% of total), with contribu- 
tions from muon neutrinos (15.58%), electron neutrinos 
(0.15%) and electron antineutrinos (0.42%). 

The neutrino flux predictions at the SciBooNE detec- 
tor location as a function of neutrino energy for both neu- 
trino and antineutrino mode running are shown in Fig.[T] 
For the low-energy BNB, the neutrino spectrum at Sci- 
BooNE and MiniBooNE are very similar, except for the 
flux normalization difference [lo| . in particular, according 
to simulations, the mean energy of the flux is expected 
to be 0.76 (0.79) GeV at the SciBooNE (MiniBooNE) de- 
tector location, reflecting the very similar fraction of i/^s 
from 7r+ and decay at the two locations. 

The systematic uncertainties in the neutrino flux pre- 
diction come from several sources in the simulation steps 
above: proton deli very /optics, secondary particle pro- 
ductions, hadronic interactions in the target or horn and 
horn magnetic field. Fig. [5] shows the flux in neutrino 
mode from tt"*" and K'^ decays, and their fractional uncer- 
tainties. Fig. [3] shows the and flux in antineutrino 
mode from tt", -k^ , and decays, and their fractional 
uncertainties. The uncertainty in tt"^ ('"'") production is 
determined from uncertainties associated with spline fits 
to the HARP 7r+ (tt") double differential cross section 
data [2j] and differences between the SW and spline fit 
central values. The HARP data used were those from 
a thin (5% interaction length) beryllium target run [8|. 
While the HARP data provide a valuable constraint on 
the BNB flux prediction, additional uncertainties result- 
ing from thick target effects (secondary re-scattering of 
protons and pions) are included through the BNB flux 
simulation. The resulting tt"*" production uncertainty is 
« 5% at the peak of the flux distribution and increases 
significantly at high and low neutrino energies. The re- 
sulting 7r~ production uncertainty is « 10% at the peak 
of the flux distribution and also increases at high and low 
neutrino energies. 

The flux from K'^ decay is dominant for > 2.0 GeV. 
Since no published data exist for production at the 
BNB primary proton beam energy, we employ the Feyn- 
man scaling hypothesis to relate production measure- 
ments at different proton beam energies to the expected 
production at the BNB proton beam energy 0] . The er- 
rors of the Feynman scaling parameters obtained from 
these measurements are then included as systematic un- 
certainties. 

Other major contributions to the flux error include un- 
certainties in the hadron interactions at the target and 
simulation of the the horn magnetic field, which both 
contribute to shape and normalization uncertainties. An 
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FIG. 1. Neutrino flux predictions at the SciBooNE detector 
as a function of neutrino energy E^, normalized per unit area, 
proton on target (POT), and neutrino energy bin width, in 
neutrino (top) and antineutrino (bottom) modes. The spectra 
is averaged within a circle with radius 2.12m from beam center 
(coincides with center of detector cross-sectional area), which 
covers the entire 3mx3m cross-sectional area of the SciBooNE 
detector. The total flux and contributions from individual 
neutrino flavors are shown. 



overall normalization uncertainty is included on the num- 
ber of protons on target (POT). All flux errors are mod- 
eled through variations in the simulation and result in a 
total error of « 7% at the peak of the flux. Quantita- 
tive constraints of each uncertainty have been determined 



FIG. 2. (Top) flux prediction at the SciBooNE detec- 
tor as a function of neutrino energy E^, in neutrino running 
mode. The total flux and contributions from tt^ and de- 
cays are shown. (Bottom) Fractional uncertainties of the z^,, 
flux prediction due to 7r+ and production from the p-Be 
interaction in neutrino running mode. The figures are from 
Ref. 
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FIG. 3. (Top) Ufj, and z?^ flux predictions at the SciBooNE 
detector as a function of neutrino energy Ei, in antineutrino 
running mode. The total flux and contributions from t/^ from 
7r~ decay and from Tr"*" and decays are shown. (Bot- 
tom) Fractional uncertainties of the i/^ and i-^, flux predictions 
due to 7r~, n'^ and K'^ production from p-Be interactions in 
antineutrino running mode with respect to the total Ufj, + 
flux. 
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from previous MiniBooNE studies 0). 



B. SciBooNE Detector 

The SciBooNE detector was located 100 m downstream 
from the berylhum target on the axis of the beam, as 
shown in Fig.lD The detector was comprised of three sub- 
detectors: a fully active and finely segmented scintillator 
tracker (SciBar), an electromagnetic calorimeter (EC), 
and a muon range detector (MRD). SciBar served as the 
primary neutrino target for this analysis. 




lOU m 44U m 



FIG. 4. Schematic overview of the Booster Neutrino Beam- 
line and the location of the SciBooNE and the MiniBooNE 
detectors. 

SciBooNE uses a right-handed Cartesian coordinate 
system in which the z axis is the beam direction and 
the y axis is the vertical upward direction. The origin 
is located on the most upstream surface of SciBar in the 
z dimension, and at the center of the SciBar scintillator 
plane in the x and y dimensions. Since each sub-detector 
is read out both vertically and horizontally, two views 
are defined: top (x vs. z projection) and side {y vs. z 
projection). 

The SciBar detector ^] was positioned upstream of 
the other sub-detectors. It consists of 14,336 extruded 
plastic scintillator strips. Each strip has a dimension of 
2.5 X 300 X 1.3 cm^. The scintillators are arranged ver- 
tically and horizontally to construct a 3 x 3 x 1.7 m"^ 
volume with a total mass of 15 tons. The dominant com- 
ponent of the SciBar detector is polystyrene (CgHg). The 
uncertainty of the total detector mass is estimated to be 
1%, including the effect of epoxy resin used to glue the 
strips. 

Each strip was read out by a 64-channel multi-anode 
photo-multiplier (MA-PMT) via a wavelength shifting 
(WLS) fiber. Charge information was recorded for each 
channel, while timing information was recorded in groups 
of 32 channels by taking the logical OR with multi-hit 
TDC modules |l3|. The timing resolution for minimum- 
ionizing particles was evaluated with cosmic ray data to 
be 1.6 ns. The average light yield for minimum- ionizing 
particles is approximately 20 photo-electrons per 1.3 cm 
path length, and the typical pedestal width is below 0.3 
photo-electron. The hit finding efficiency evaluated with 
cosmic ray data is more than 99.8%. The minimum 
length of a reconstructable track is approximately 8 cm 
(three layers hit in each view). The track finding effi- 
ciency for single tracks of 10 cm or longer is more than 



99%. Large samples of events were visually scanned to 
ensure that the track finding algorithm is efficient, has 
high purity, and is unbiased. The scanning showed that 
the purity of the track-finding algorithm using the in- 
formation available from the detector was efficient for 
separating events into SciBar 1, 2, 3-Track samples. 

The EC is located just downstream of SciBar, and is 
designed to measure the electron neutrino contamination 
in the beam and tag photons from 7r° decay. The EC is 
a "spaghetti" type calorimeter comprised of 1 mm diam- 
eter scintillating fibers embedded in lead foil [l3|. The 
calorimeter is made of 64 modules of dimensions 262 x 
8x4 cm"^. The fibers are bundled in two independent 
groups of 4 X 4 cm^ transverse cross section, read at both 
ends by Hamamatsu PMTs. 

The MRD was installed downstream of the EC and is 
designed to measure the momentum of muons produced 
by CC neutrino interactions. It comprised of 12 iron 
plates with thickness 5 cm sandwiched between planes 
of 6 mm thick scintillation counters; there were 13 al- 
ternating horizontal and vertical planes read out via 362 
individual 2 inch PMTs. Each iron plate measured 274 x 
305 cm^. The MRD measured the momentum of muons 
up to 1.2 GeV/c using the observed muon range. Charge 
and timing information from each PMT were recorded. 
The average hit finding efficiency is 99%. 



C. Detector Response Simulation 

The GEANT4 framework is used for the detector simu- 
lation. The Bertini cascade model within GEANT4 [ll] 
is used to simulate the interactions of hadronic parti- 
cles with detector materials. The detector simulation 
includes a detailed geometric model of the detector, in- 
cluding the detector frame, experimental hall, and soil. 
The energy loss of a particle in each single SciBar strip 
and each individual EC sensitive fiber is simulated. The 
energy deposition is converted in the detector response 
taking into account the Birk's saturation of the scintil- 
lator, the light attenuation along the fibers, the Poisson 
fiuctuation of the number of photoelectrons, the PMT 
resolution, and electronic noise. The crosstalk in nearby 
SciBar channels is also simulated. 

In SciBar the timing of each hit is simulated from the 
true time of the corresponding energy deposition, cor- 
rected by the travel time of the light in the WLS fiber 
and smeared by the timing resolution. 

For the detector simulation of the MRD, true energy 
deposition in each scintillator is converted to ADC counts 
using the conversion factor measured with cosmic muons. 
The attenuation of light in the scintillator as well as elec- 
tronics noise are simulated. The time of energy deposi- 
tion is digitized and converted into TDC counts. 

The input parameters of the detector simulation are 
derived from laboratory measurements and calibration 
data. The features of the simulation have been systemat- 
ically compared and tuned with cosmic ray and neutrino 
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data. 

A more detailed description of the detector simulation 
is given in (l6j . 

D. Neutrino Interaction Simulation 

In the SciBooNE experiment, neutrino interactions 
with carbon and hydro gen in the SciBar detector are 
simulated by the NEUT and NUANCE flQ^ pro- 

gram libraries. NEUT is used in Kamiokande, Super- 
Kamiokande, K2K, and T2K experiments, while NU- 
ANCE is used in MiniBooNE. Several Monte Carlo sam- 
ples with different NEUT and NUANCE implementa- 
tions are produced, and compared to the SciBooNE neu- 
trino data. For the analysis presented in this paper we 
use NUANCE as the neutrino interaction simulation code 
because this Monte Carlo matches what was used in the 
MiniBooNE oscillation analysis ^j. The same analysis 
will also be repeated using NEUT and presented in Ap- 
pendix |B1 The total number of protons on target (POT) 
collected in neutrino mode is 0.99 x 10'^" while the POT 
for the antineutrino mode is 1.51 x 10^°. The expected 
number of events in the SciBooNE detector for each neu- 
trino (antineutrino) interaction is listed in Tab. HI 



TABLE I. The expected number of events in each neutrino 
interaction estimated by NUANCE at the SciBooNE detec- 
tor location with the neutrino beam exposure of 0.99 x 10^" 
protons on target for neutrino mode and of 1.51 x 10^" for 
antineutrino mode. The 9.8 ton fiducial volume of the SciBar 
detector is assumed. Charged Current and Neutral Current 
interactions are abbreviated as CC and NC, respectively. 



Interaction 



NUANCE 
# Events # Events 

neutrino mode antineutrino mode 



CC QE 
CC single- TT 
CC coherent tt 
CC DIS-hOther 
NC 



45,163 
24,437 
1,706 
3,049 
29,118 



15,361 
6,413 
1,326 

I, 518 

II, 686 



elastic scattering {uN — iN'hadrons), where A'' and N' 
are the nucleons (proton or neutron), i is the lepton (elec- 
tron, muon or neutrino), and m is the meson. NUANCE 
also simulates resonantly-produced multi-pion and kaon 
production (AT+A, AT+E). In nuclei, interactions of the 
mesons and hadrons with the nuclear medium are simu- 
lated following the neutrino interactions. 

In addition to neutrino interactions inside SciBar, we 
also simulate interactions in the EC and MRD and the 
surrounding materials (the walls of the detector hall and 
soil). 



1. Quasi-elastic Scattering 

The dominant interaction in SciBooNE is CC quasi- 
elastic scattering, which is implemented using the Smith 
and Moniz model [2lj . The nucleons are treated as quasi- 
free particles with Fermi motion and Pauli blocking taken 
into account. The Fermi surface momentum (pf) for car- 
bon is set to 220 MeV/c and the nuclear potential (Eb) 
is set to 34 MeV/c, as extracted from electron scatter- 
ing data [2^ taking into account neutrino vs. electron 
scattering differences jlSl- For the vector form factor, 
NUANCE uses the BBA-2003 form factor [H with a 
dipole form for the axial form factor with an adjustable 
axial mass, M^^ for Carbon (M^^-C) and for Hydrogen 
(M^^-H). The M^^-H is used only in the antineutrino 
mode running. In NUANCE, an empirical Pauli-blocking 
parameter, k is introduced [25l | to better describe the 
MiniBooNE quasi-elastic data at low momentum trans- 
fer. When K > 1, the Pauli-blocking of final state nu- 
cleons is increased and hence the cross section at low 
momentum transfer is suppressed. The values of {AI^^- 
C = 1.234, M^^-H = 1.0) GeV/c2 and k = 1.022 are 
used [23| . The same Fermi momentum distribution and 
nuclear potential are used in all other neutrino-nucleus 
interactions except for coherent tt production. 



S. Meson Production via Baryon Resonances 



For neutrino mode beam exposure, the total number of 
CC interactions predicted by NUANCE integrated over 
the SciBooNE flux in the 9.8 ton SciBar fiducial volume 
is 7.44 X 10"* for 0.99 x lO^o POT. For antineutrino mode 
beam exposure, the total number of CC interactions pre- 
dicted by NUANCE integrated over the SciBooNE flux 
in the 9.8 ton SciBar fiducial volume is 2.46 x lO"* for 
1.51 X 10^0 POT. 

The targets handled in NUANCE are proton, neutron, 
and carbon nuclei. The types of neutrino interactions 
simulated in both Neutral Current (NC) and CC are: 
elastic and quasi-elastic scattering {vN — >■ tN'), single 
meson production {vN — >■ £N'm) , single gamma produc- 
tion {vN — £N'j), coherent tt production, and deep in- 



The second most frequent interaction in SciBooNE is 
the resonant production of single pion, kaon, and eta 
mesons which NUANCE described with the model of 
Rein and Sehgal [2^]. The Rein and Sehgal model as- 
sumes an intermediate baryon resonance, A^*, in the 
reaction of £N*,N* N'm. All intermedi- 

ate baryon resonances with mass less than 1.7 GeV/c^ 
are included. Interactions with invariant masses greater 
than 2 GeV/c^ are simulated as deep inelastic scatter- 
ing. A re-interactions (AA^ — > A^A^) which do not lead 
to a mesonic final state are also simulated. This re- 
interaction probability is assumed to be 10% (20%) for 
A++/- ( A+/0) resonances. To determine the angular 
distribution of final state pious, the method of [27| is 
used for the P33(1232) resonance. For other resonances. 
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the directional distribution of the generated pion is cho- 
sen to be isotropic in the resonance rest frame. The 
axial- vector form factors are formalized to be dipole with 
MX"" = 1.10 GeV/c2. 



3. Coherent Pion Production 

Coherent pion production is a neutrino interaction 
with a nucleus which remains intact, releasing one pion 
with the same charge as the incoming weak current. Be- 
cause of the small momentum transfer to the target nu- 
cleus, the outgoing pion tends to be emitted in the for- 
ward direction, closely following the incoming neutrino 
direction. The formalism developed by Rein and Seh- 
gal 0, [2^ is used to simulate such interactions. The ax- 
ial vector mass, M^°'''='''="* , is set to 1.03 ± 0.28 GeV/c^. 
The total and inelastic pion-nucleon cross sections from 
the original Rein-Sehgal publications [1^ are ob- 
tained from fits to PDG data and implemented as a func- 
tion of pion energy. The Rein and Sehgal model predicts 
the CC coherent 7r+ production rate to be approximately 
1% of the total neutrino CC rate in SciBooNE. 



4- Deep Inelastic Scattering 

The deep inelastic scattering (DIS) cross section is 
calculated using the GRV98 parton distribution func- 
tions [131 • Additionally, we have included the corrections 
in the small region developed by Bodek and Yang [3l| . 
The DIS contribution slowly increases for W values start- 
ing at 1.7 GeV and becomes the only source of neutrino 
interactions above W>2 GeV. This is done to create a 
smooth transition between the resonance and DIS mod- 
els and ensure continuity in distributions of kinematics 
and hadron multiplicity in the region of overlap. Tab. |TT] 
summarizes the parameter choices used in NUANCE. 



TABLE II. NUANCE parameters used for neutrino interac- 
tion simulation. 



PF 220 MeV 

Eb 34 MeV 

M^^-C 1.234 GeV 

M^^-H 1.0 GeV 

K 1.022 

MX" 1.10 GeV 

j^coherent ^ Qg Q^y 

M^"" 1.3 GeV 



5. Intra-nuclear Interactions 

Following production, the intra-nuclear interactions of 
mesons and nucleons are simulated using a cascade model 
in which the particles are traced until they escape the 
target nucleus. 

Although we only use kinematic information from the 
final state muon in this analysis, the simulation of inter- 
nuclear interaction is important since the pions/protons 
emitted from the nucleus can affect the selection criteria 
and bias the kinematics of the final state muon. 

Inelastic scattering, charge exchange and absorption of 
pious in nuclei are simulated. For inelastic scattering and 
charge exchange interactions, the direction and momen- 
tum of pions are affected. In the scattering amplitude, 
Pauli blocking is also taken into account. A more detailed 
description of the intra-nuclear interaction simulations in 
NUANCE can be found in [ll|. 

III. NEUTRINO MODE ANALYSIS 

In this analysis, we use a sample of high energy 
events identified by muons penetrating the SciBar, EC, 
and MRD detectors. According to our Monte Carlo sim- 
ulation, 43% of these high energy neutrinos come from 
K'^ decays in the beam. 

The entire neutrino-mode data set for the SciBooNE 
experiment is used for this analysis. The neutrino run 
occurred from October 2007 to April 2008. A total of 
0.99 X lO'^" POT, after all data quality cuts, was collected 
with an efficiency of 95%. 

The hit threshold in SciBar is set at 2.5 photoelec- 
trons (corresponding to roughly 0.25 MeV) for a scintil- 
lator strip. The hits in each view are then associated into 
two-dimensional (2D) tracks in each view using a cellu- 
lar automaton algorithm developed in K2K [S^]- First, 
a correction for cross talk is applied to both data and 
MC prediction before 2D track reconstruction. Then, 
the remaining hits are categorized into "clusters" , where 
a cluster is one or more hits in adjacent scintillator strips. 
Segments are formed, which connect individual clusters 
no more than one scintillator strip apart. Any segments 
which share a cluster are then connected, provided the 
of a least squares linear fit remains acceptable. Three 
dimensional tracks are formed from 2D tracks in both 
views by requiring the timing between the 2D tracks to 
be within 50 ns, and the start and end point in the z 
direction for the 2D tracks to be within 6.6 cm. 

A reconstructed muon track (referred to as a SciBar- 
MRD matched track) consists of a SciBar reconstructed 
track matching a reconstructed track in the MRD. The 
fiducial volume (FV) is applied requiring the upstream 
edge of the track to be within |a;| < 130 cm, \y\ < 130 
cm, and 5.24 < < 149.34 cm. This corresponds to 
a 9.7 m^ FV and 9.8 tons of fiducial mass. This selec- 
tion criterion is not the same as is used in previous Sci- 
BooNE pubhcations of 10.6 tons [Ullllll]; the selection 
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has been optimized for the particular analysis described 
here. The difference is due to a more conservative cut in 
the z direction to better reduce background events and 
more accurately reconstruct kinematic variables associ- 
ated with the SciBar tracks. 

As only MRD events are selected, a muon-hypothesis 
track selection is unnecessary, and the FV selection has 
been modified to boost the number of events in the sam- 
ple without adding backgrounds. The track is also re- 
quired to be in time with the neutrino beam, i.e. within 
< t < 2, 000 ns. The extrapolation of the SciBar track 
is required to be matched with a track in MRD within 30 
cm and within a 100 ns time difference. The matching 
requires a track that penetrates at least four consecutive 
MRD layers. If no MRD track is found, MRD hits can 
be used for the matching. The MRD hits, which in this 
case are required to be within a cone with an aperture of 
±0.5 rad, should have a time within 100 ns with respect 
to the SciBar track and should be within 10 cm of the 
extrapolated SciBar track on the MRD first layer. 

The SciBar-MRD matched track can be further classi- 
fied based on the end point of the MRD matched track: 
if the SciBar-MRD matched track stops in the MRD, the 
track is defined as SciBar-MRD stopped. If the track goes 
through all of the MRD planes, the track is classified as 
SciBar-MRD penetrated. If the SciBar-MRD matched 
track escapes from the side of the MRD before the most 
downstream plane, the event is classified as SciBar-MRD 
side-escaped. 

In this analysis we will consider only events with one 
SciBar-MRD penetrated track. This defines a sample en- 
riched in high energy muon neutrino events. The selected 
muons will have a minimum momentum of 1.0 GeV/c and 
muon angle relative to beam axis of less than 45 degrees. 

The data sets for SciBooNE are subject to contami- 
nation by cosmic backgrounds (due to lack of detector 
shielding or veto) . The vast majority of the cosmic back- 
grounds come from cosmic muons, which can mimic a 
muon signal from a interaction in SciBar. The selec- 
tion cuts remove most of the cosmic background muons 
due to differences between cosmic background muons 
and muons created by neutrino beam interactions. The 
FV and timing cuts in SciBar remove the cosmic back- 
grounds coming from outside the detector and beam win- 
dow. Most of the cosmic backgrounds enter the detec- 
tor from above, producing vertical tracks (neutrino beam 
event tracks are generally horizontal). The vertical cos- 
mic tracks do not pass through both SciBar and MRD 
and no SciBar-MRD matched track can be reconstructed, 
failing the cut requiring a SciBar-MRD matched track. 
For events with a SciBar-MRD matched track, the cos- 
mic background contamination in the beam timing win- 
dow is 0.5%, estimated using a beam-off timing window 
(5,000 < t < 15,000 ns). To further correct for the 
cosmic backgrounds, identical selection cuts are made in 
data in the timing region 5, 000 < t < 15, 000 ns, far 
outside the beam timing window and 5 times as long for 
better statistics. The data numbers and distributions se- 



lected in this cosmic background region are then divided 
by 5 and subtracted from the corresponding data events 
and distributions selected in the beam timing window. 
The resulting cosmic background subtracted data num- 
bers and distributions assume no further cosmic back- 
ground contamination. There is no cosmic background 
modeling in the MC. The number of events from cos- 
mic backgrounds, neutrino interactions in the EC and 
MRD, and events with neutrino interactions in the ma- 
terial surrounding the SciBar detector are negligible. For 
all three samples combined, the number of background 
events is estimated to be 26 events: 15 from cosmic, 9 
from backscattering neutrino events in EC/MRD, and 2 
from interaction outside the SciBooNE detectors. 

The high energy muons penetrate through the entire 
SciBooNE detector so the reconstruction of the total 
muon energy could not be done. The reconstructed muon 
angle relative to beam axis will be used as the primary 
kinematic variable for the analysis. Neutrinos produced 
from decay have a higher energy on average than 
neutrinos from 7r+ decay. Therefore, the angular distri- 
bution of the resulting muon from the neutrino interac- 
tion of a neutrino from will be more forward peaked 
than those from neutrinos from 7r+. 

After all the previously mentioned selection cuts, the 
resulting data and MC samples are further divided into 
three separate samples based on the number of SciBar re- 
constructed tracks in each event. The number of SciBar 
reconstructed tracks for data and MC (along with its dif- 
ferent contributions) can be seen in Fig. [51 The SciBar 1- 
Track sample contains mostly (78%) CCQE interactions, 
the SciBar 2-Track sample contains both CCQE(46%) 
and CCItt (42%) and the SciBar 3-Track sample contains 
mostly CCl7r(62%) events together with QE and DIS in- 
teractions. The reconstructed muon angle distributions 
for the SciBar 1, 2 and 3-Track samples are shown in 
Fig. ini Fig. ini shows an excess in the MC events with re- 
spect to the data for low angle muon reconstructed tracks 
which is largely attributed to the difficulty in correctly 
modeling nuclear final state interactions. We believe that 
this MC / data discrepancy at low reconstructed muon an- 
gle is due to modeling of proton emission; the effect is 
covered by the uncertainty included for the modeling in 
Sec. El 

Tab, mil shows the number of events coming from , 
TT+, data, and MC (NUANCE) along with efficiencies 
and purities (defined in Eqs. |2l (3]) for the three selected 
samples. 

We define the efficiency and purity using MC simulated 
events as: 

# selected events 

# generated ' 

# true in selected sample 

# selected events ' 

The mean energy and mean angle (with respect to pro- 
ton beam direction) of the selected and the mean 
energy for from the selected in each of the three 



9 



-^Data 
— MC 



Number of Reconstructed SciBar Tracks 



FIG. 5. Number of SciBar reconstructed tracks for the se- 
lected event sample in neutrino mode for data and NUANCE 
MC. MC histogram includes all events and shows the 
component of from . 



TABLE III. Number of events for data and for NUANCE MC 
after selection cuts for the neutrino mode analysis. MRD in- 
dicates the SciBar-MRD sample, and MRD pen. stands for 
SciBar-MRD penetrated sample. The cosmic backgrounds 
have been subtracted from the data. The first column con- 
tains the number of predicted events with a coming from 
while the second column contains the number of pre- 
dicted events with a coming from tt"*". The third column 
and fourth column represent the number of events in data 
and MC and the last two columns are the efficiency {e[K^)) 
and purity {■k{K^)) for i/^ from in MC. The total event 
prediction from simulation is labeled as MC. 



Event Sel. 






Data 


MC 


<K+) 


■k(K+) 


MRD 


2,867 


18,173 


27,049 


22,142 


58% 


13% 


MRD pen. 


1,508 


1,700 


3,365 


3,286 


31% 


46% 


Single 


1,313 


1,666 


3,188 


3,053 


27% 


43% 


1-Trk 


509 


1,159 


2,050 


1,723 


10% 


30% 


2-Trk 


497 


438 


834 


950 


10% 


52% 


3-Trk 


189 


57 


206 


250 


4% 


76% 



samples is summarized in Tab. IIVI Fig. [7] shows the 2- 
dimensional distribution of if + production angle relative 
to proton beam axis versus true production energy 
for predicted K'^ events selected in the SciBar 1-Track, 
2-Track and 3-Track MC samples for the neutrino mode 
analysis. 

The background from K~ and K'^ decays are very 
small. The expected number of and events pre- 
dicted for all three samples combined are 5 and 3 events, 
respectively. 
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(a) 1-Track Sample 
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(b) 2-Track Sample 




Reconstructed Muon Angle (degrees) 

(c) 3-Track Sample 

FIG. 6. Reconstructed muon angle for the SciBar 1-Track, 
2-Track, and 3-Track samples in the neutrino mode analy- 
sis. The background contributions from K~ and are very 
small but included in the MC histogram. The grey area rep- 
resents the total systematic uncertainty in the MC. 



TABLE IV. NUANCE MC predicted mean energy and mean 
angle (with respect to the proton beam direction) for the se- 
lected samples and predicted MC energy for the selected 
i/fi from for the three SciBar samples in neutrino mode 
running. Errors correspond to the error on the mean values. 
RMS of the relative distributions is reported between squared 
parenthesis. 



IV. ANTINEUTRINO MODE ANALYSIS 

The dataset used in this analysis consists of 
antineutrino-mode events collected between June 2007 
and August 2008 corresponding to a total of 1.51 x 10^° 
POT, after all data quality cuts. The data collection 



Ek+ [RMS](GeV) | 6>y+ [RMS](deg) | j;,^[RMS](GeV) 

1- Trk 3.6±0.1[1.2] 4.3±0.1[2.1] 3.0±0.1[1.0] 

2- Trk 378^0. 1[1. 2] " 4.1±0.1[2.0] 3.3±0.1[1.0] 
FTrk OzfcO.ljl.l] 3.9±0.1[1.9^ 3.5±0.1[0.9] 
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(a) 1-Track Sample 
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(c) 3- Track Sample 

FIG. 7. Two-dimensional plots of the energy vs angle relative 
to the proton beam axis for the selected events in the 
SciBar 1-Track, 2- Track, and 3- Track MC samples for the 
neutrino mode analysis. 



efficiency was also 95% in this case. 

The majority of events in antineutrino mode running 
are from 7r~ decay, with a significant i/^ contribution 



coming from the decay of positively charged tt"*" and . 
At higher energies (above 2 GeV), the positively charged 
"wrong-sign" particles, mainly and ^ are strongly 
boosted in the forward direction, traveling downstream 
nearly parallel to the beam axis direction, such that they 
are not defocused by the magnetic horn. These particles 
create a 40% background of wrong sign neutrino inter- 
actions in the antineutrino mode run. In this analysis, 
we select a sample of high energy i/^ coming from . 
The analysis strategy is very similar to what is described 
in Sec. IIIII with the difference that the background com- 
position in the selected sample is different. Background 
events result from high energy from tt"*" produced in 
the initial p-Be interaction, and the decay of 7r~ (also 
produced in the initial p-Be interaction) that generate 
high energy D^. Both the and result in high en- 
ergy muons (with negative and positive charge respec- 
tively) that penetrate the SciBar, EC and MRD detec- 
tors. Positive and negative muons are indistinguishable 
in SciBooNE due to the lack of magnetic field in our de- 
tector system. 

As in the neutrino mode analysis, events passing the 
base selection cuts are further divided into three sam- 
ples based on whether the events contain 1, 2 or 3 
SciBar reconstructed tracks. The number of SciBar re- 
constructed tracks for data and MC (along with the par- 
ent) can be seen in Fig. [51 The SciBar 1-Track sample 
contains mostly (81%) charged current quasi-elastic in- 
teraction events. The SciBar 2-Track sample is evenly 
split between charged current quasi-elastic (40%) and 
charged current resonant pion interactions (44%), with 
a tiny contribution from charged current multi-7r/DlS 
interactions (6%). The SciBar 3- Track sample contains 
mostly charged current resonant pion interactions (61%), 
with small contributions from charged current quasi- 
elastic (12%) and charged current multi-7r/DIS interac- 
tions (21%). The reconstructed muon angle distributions 
for the SciBar 1, 2 and 3- Track samples are shown in 
Fig. [9l In this case the events coming from have 
high energies, peaking at smaller angles, while the lower 
energy 7r+ and 7r~ background distributions are spread 
more evenly across a broader range of angles. 

In this analysis, we also found a slight disagreement 
between data and Monte Carlo at low angle regions for 
the SciBar 2 and 3-Track samples for the reconstructed 
muon angle (see Fig. [9|). This disagreement is caused, as 
in the analysis, by the imperfect modeling of nuclear 
and intra-nuclear interactions in our neutrino interaction 
simulation program. 

Tab. |V] summarizes the selected number of events in 
Data and MC (NUANCE) for the three selected samples, 
together with the efficiency and purity of (defined in 
Eqs.iaEl) . 

The mean energy and mean angle (with respect to the 
proton beam direction) for the selected K'^ and the mean 
energy for from the selected in each of the three 
samples is summarized in Tab. IVII Fig. [10] shows the 2- 
dimensional distribution of if + production angle relative 
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FIG. 8. Number of SciBar reconstructed tracks for the se- 
lected event sample in antineutrino mode for data and NU- 
ANCE MC. MC histogram includes all events and shows 
the component of v^, from and + vr"*" shows the com- 
ponent of Vfj, from and Tr"*" . 



TABLE V. Number of events for data and for NUANCE MC 
after selection cuts for the antineutrino mode analysis. MRD 
indicates the SciBar-MRD sample, and MRD pen. stands for 
SciBar-MRD penetrated sample. The cosmic backgrounds 
have been subtracted from the data. The first column con- 
tains the number of predicted events with Ufj, coming from , 
the second column contains the number of predicted events 
with coming from tt"*", and the third column contains the 
number of predicted events with i/^ coming from ■k~ . Fourth 
column and fifth column represent the number of events in 
data and MC and the last two columns are the efficiency 
{e{K+)) and purity {tv{K+)) for ly^ from K+ in MC. The 
total event prediction from simulation is labeled as MC. 



Event Sel. 








Data 


MC 






MRD 


705 


2,287 


6,167 


11,528 


9,499 


52% 


7% 


MRD pen. 


326 


385 


698 


1,790 


1,452 


24% 


22% 


Single fi 


283 


375 


691 


1,728 


1,389 


21% 


20% 


1-Trk 


119 


230 


589 


1328 


965 


9% 


12% 


2-Trk 


103 


121 


83 


296 


317 


8% 


32% 


3-Trk 


36 


19 


15 


75 


74 


3% 


49% 
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FIG. 9. Reconstructed muon angle for the SciBar 1-Track, 2- 
Track, and 3- Track sample in the antineutrino mode analysis. 
The background contributions from K~ and are small but 
included in the MC histogram. The grey area represents the 
total systematic uncertainty in the MC. 



to proton beam axis versus true production energy 
for predicted events selected in the SciBar 1-Track, 2- 
Track and 3-Track MC samples for the antineutrino mode 
analysis. As can be seen in Fig. [TOl the are much 
more forward than in neutrino mode (Fig. [7|) otherwise 
they would have been swept out by the horn. 

The number of neutrinos coming from K~ and 
decays are small: the number of expected and if^ 
events predicted for all three samples combined are 25 
and 5 events, respectively. Fig. [TT] shows the recon- 
structed muon angle for the 2- Track sample with the 
and Kl contributions. Interactions from cosmic back- 
grounds, neutrino interactions in the EC and MRD, and 
events with neutrino interactions in the material sur- 
rounding the SciBar detector make a negligible contribu- 
tion in this analysis as well. For all three samples com- 



TABLE VI. NUANCE MC predicted mean energy and mean 
angle (with respect to the proton beam direction) for the se- 
lected K'^ samples and predicted MC energy for the selected 
i/fi from for the three SciBar samples in antineutrino mode 
running. Errors are the error on the mean values. RMS of the 
relative distributions is reported between squared parenthesis. 





Ei^.+ [RMS](GeV) 


ej^+[RMSl(dcg) 


i;4RMS](GeV) 


1-Trk 


4.1±0.1[1.2] 


2.4±0.2[1.8] 


3.1±0.1[1.1 




2-Trk 


4.4±0.1[1.2] 


1.9±0.1[1.4] 


3.4±0.1[1.0 




3-Trk 


4.6±0.2[1.1] 


1.6±0.2[1.0] 


3.6±0.2[1.0 





bined, 19 events are estimated from cosmic backgrounds. 
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FIG. 10. Two-dimensional plots of the energy vs angle relative 
to the proton beam axis for the selected events in the 
SciBar 1-Track, 2- Track, and 3- Track MC samples for the 
antineutrino mode analysis. The are much more forward 
than in neutrino mode (Fig. [7| otherwise they would have 
been swept out by the horn, this can be seen also in the 
angular distribution difference between Tab. II VI and Tab. I VII 



1 event is estimated from backscattering neutrino events 
in EC and MRD, and 3 events are estimated from neu- 
trino interactions outside the SciBooNE detector. 




Reconstructed Muon Angle (degrees) 



FIG. 11. The various meson contributions that decay into 
v,j, and Uf_i as a function of reconstructed muon angle in the 
2-Track sample in antineutrino mode. 



V. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are included by using a co- 
variance matrix that includes the correlated and uncor- 
related errors among the muon angle bins in the SciBar 
1, 2, 3- Track samples. The covariance matrix is deter- 
mined by calculating the correlated event changes for a 
given systematic uncertainty and then combining these 
(the different systematic uncertainties are assumed to be 
uncorrelated) to form the matrix. Statistical errors are 
also included as uncorrelated terms in the diagonal ele- 
ments. 



A. Neutrino Beam Uncertainties 

The uncertainties in the delivery of the primary pro- 
ton beam to the beryllium target and the primary beam 
optics, secondary hadron production in proton-beryllium 
interactions, hadronic interactions in the target and horn 
and the horn magnetic field model are considered. The 
change in the neutrino beam spectrum due to these un- 
certainties is calculated by drawing random parameter 
vectors (all the beam systematics are varied within their 
uncertainties when drawing random parameter vectors) 
and weighting each event by a factor corresponding to 
the variation of the yield of the parent meson with the 
given momentum and angle. 

To evaluate the neutrino beam uncertainties for both 
the neutrino and antineutrino mode running, a thousand 
random parameter vectors are generated, resulting in a 
thousand neutrino beam flux predictions and a thousand 
individual weights for each MC event. The thousand in- 
dividual sets of weights for the MC simulation are passed 
through each analysis resulting in a thousand individual 
outcomes for each analysis. The correlated event changes 
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associated with the thousand individual outcomes are 
used to form the error matrix for the beam systematic 
uncertainties. 

For K~ production error, which is not accounted for in 
the previous neutrino beam uncertainty weights, a con- 
servative 100% production uncertainty is apphed. 

B. Detector Uncertainties 

1. PMT Crosstalk and Resolution 

The crosstalk of the MA-PMT was measured to be 
3.15% for adjacent channels, with an absolute error of 
0.4% [l^. The single photoelectron resolution of the 
MA-PMT is set to 50% in the simulation, and the ab- 
solute error is estimated to be ±20%. Several complete 
MC simulation sets, with the crosstalk level and single 
photoelectron resolution separately varied within their 
uncertainties, are prepared and the changes in the final 
results using the varied MC simulation sets are taken as 
the systematic uncertainties. 

2. Scintillator Quenching 

Birk's constant for the SciBar scintillator was mea- 
sured to be 0.0208 ± 0.0023 cm/MeV [16.] and is varied 
within the measurement to evaluate the systematic error. 

3. Hit Threshold 

The conversion factors from the ADC counts to the 
photoelectron were measured for all 14,336 MA-PMT 
channels in SciBar. The measurement uncertainty was 
at the 20% level. Since the number of photoelectrons 
(2.5 p.e.) for the SciBar hit threshold was used, the hit 
threshold is varied by ±20% to evaluate the systematic 
error for SciBar track reconstruction. 



4. TDC Dead time 

The TDC dead time is set to 55 ns in the MC simula- 
tion, with the error estimated to be ±20 ns [s^. A MC 
simulation set with variations in the TDC dead time was 
prepared to determine the systematic error. 

C. Cross section and nuclear model uncertainties 

1. CC Quasi Elastic Scattering Cross Section 

A systematic uncertainty of ±0.234 GeV is assigned to 
the M^^-C to span the difference between the value used 
and the value from the global fit to historical data [35| . 



The difference between k = 1.000 and k = 1.022 is 
also assigned as a systematic uncertainty and added in 
quadrature to the M^^-C error. M^^-H error gives a 
negligible contribution. 

2. CC Resonant Pion Production Cross Section 

From a previous K2K measurement (36| . the uncer- 
tainty in the resonant pion scattering cross section is esti- 
mated to be ±20%. An additional uncertainty is assigned 
to account for the observed disagreement between the 
SciBooNE CC l7r-enriched data samples and the MC fll| . 
The uncertainty is evaluated by re- weighting CC reso- 
nant pion events as a function of true such that they 
match the observed distribution in the SciBooNE data. 
Resonance decays leading to multi-pion final states are 
also included in the model and are simulated assuming 
M^'^ = 1.30 GeV/c^. This value of M^'^ is chosen 
strictly to ensure that the total CC cross section pre- 
diction reproduces previous experimental data. 

3. Multi Pion Production and Deep Inelastic Scattering 
Cross Section 

The uncertainties for multi pion and deep inelastic 
scattering (DIS) cross sections are set respectively at 
±40%. The 40% multi-pion and DIS uncertainties come 
from a comparison of the NUANCE predictions to the 
existing multi-pion data in [33, [11] . 

4- Pion Interaction in the Initial Target Nucleus 

For pions produced by neutrino interactions, uncer- 
tainties on the cross sections for pion absorption, pion 
inelastic scattering, and pion charge exchange in the nu- 
cleus are considered. The values of absorption, inelastic 
scattering, and charge exchange for pions are varied by 
±30% in the MC independently while keeping the total 
number of neutrino events with pion processes fixed to 
determine the systematic uncertainties of each process. 

5. Proton Emission 

A systematic uncertainty associated with the emission 
of a proton following pion absorption in the nucleus is 
also included. The proton emission following a pion ab- 
sorption alters both the number of reconstructed tracks 
and the event kinematics. Turning off this nuclear pro- 
cess in the MC results in better agreement between the 
data and the MC, thus, indicating that the modeling may 
be inaccurate. We have estimated the uncertainty asso- 
ciated with proton emission by calculating the difference 
for the MC with this process turned on and off. The 
uncertainty is calculated for both NUANCE and NEUT, 
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with NEUT giving the larger systematic uncertainty, so 
the uncertainty is conservatively estimated using NEUT 
for both MCs. 



VI. FITS TO DETERMINE THE K+ 
PRODUCTION AND RATE NORMALIZATION 

The final state reconstructed track multiplicity and 
the muon angular distributions for the selected high en- 
ergy events are different for neutrinos produced by pion 
and kaon decays as presented in Sec. IIIII and in Sec. IIVI 
Therefore, we fit these distributions in order to isolate 
neutrinos from kaon decays and minimize the following 
function: 

2 _ 2 I 2 _ 

% Xu ' Xu 

N 

i,j 
M 

p,1 

Cross section values (CCQE and CCItt for v and v) 
and uncertainties are included in the fit as described re- 
spectively in Appendix VK\ and in Sec. IVl 

Two separate normalizations with respect to the 
beam MC predictions (Eq. [1} are extracted from the 
selected samples, a K^^od^ production flux at the 

beam target, and a K^ate^ the K^^od^ neutrino cross- 
sections. These two types of determinations are needed 
since the K^^^^ normalization can be used to determine 

a measured production cross section and the if^^e 
normalization can be used by BNB experiments as better 
estimates of the production rate in their beam. 

The function in Eq. |4] contains two terms: the for- 
mer xt term is associated with events for neutrino mode 
running and the latter Xp term is associated with events 
for antineutrino mode running. The two functions 
are assumed to be effectively uncorrelated since the cross 
section uncertainties for the antineutrino mode data are 
small compared to the statistical and background uncer- 
tainties. In the neutrino mode running, all three samples 
(SciBar 1, 2, 3-Track) are used simultaneously in the fit 
including their correlated bin-to-bin uncertainties. Only 
bins with 10 or more events are included in the x^. N 
and M are the number of bins used in the three recon- 
structed angle distributions in neutrino and antineutrino 
mode, respectively. -^1°^) and are the numbers of 

observed and predicted events in the i(j)-th angle bin 
for the neutrino mode analysis. Af and M^J'^s are the 

same quantities for the antineutrino in the p(g)-th angle 
bin. 

For the if + production analysis, the functions that de- 
scribe the number of predicted events N^^^J^ and M^T^!'' 



are given by Eq. lAll for neutrinos and Eq. IA3I for an- 
tineutrinos in Appendix Rl {V^yg)ij and {V^yg)pq are the 
elements of the covariance matrix for neutrino and an- 
tineutrino mode for each of the systematic uncertainties 
described in Sec. El V^^^f (Ktot) represents the statistical 
error in neutrino mode running (antineutrino mode run- 
ning). An initial neutrino mode and antineutrino mode 
combined x^ minimization is performed to determine the 
best cross-section normalization factors for both i/^ and 
as described in Appendix \X\ Pull terms on the cross- 
sections normalization factors are added to keep the min- 
imization physical. After the initial combined x^ min- 
imization, the cross-section weights are fixed in Nf^'^'^ 
and M^^^'^ at their minimized values and the pull terms 
are removed from the x^ to evaluate the total (statisti- 
cal-|-systematical) uncertainty on the production or 
rate. These cross-section weights are initially minimized 
to allow for better agreement between data and MC in 
the plots and do not affect either K'^ production or rate 
weights because the large uncertainties on the neutrino 
cross-section values are already taken into account in the 
covariance matrix. 

TABLE VII. Summary of the fit results for the cross-section 
normalization factors as described in Appendix[X]with respect 
to NUANCE predictions. The top four values are cross section 
normalization values for and coming from tt^ and K'^ 
while the last two are normalization factors for coming 
from 7r~. 



TT+/K + 




Fit Value 




CCQE 


in 


1/ 


mode 


1.17±0.14 




CCQE 


in 


V 


mode 


1.07±0.25 




CCItt 


in 


V 


mode 


0.89±0.25 




CCItt 


in 


v 


mode 


0.91±0.26 


TV~ 








CCQE 


in 


V 


mode 


1.50±0.21 




CCItt 


in 


V 


mode 


1.49±0.29 



The rate is measured by minimizing the same x^ 
function as described in Eq. |3] but using Eq. IA2I and 
Eq. 1^4] respectively for neutrino and antineutrino to pre- 
dict the number of events. The covariance matrices used 
for the rate measurement does not include the cross 
section uncertainties for the from in contrast to 
the matrices used for the flux measurement. A sum- 
mary of the cross section normalizations is presented 
in Tab. IVIII Many of these values are consistent with 
low energy precision cross-section measurements from the 
MiniBooNE experiment 0, l39l - l4l| though the two sets 
of cross-section values are measured at different energies. 
The MiniBooNE collaboration measures cross-sections at 
neutrino energies less than 2 GeV while the cross-section 
values listed in Tab. IVlJ are at neutrino energies greater 
than 3 GeV (as could be seen from Tab.|lV]and Tab.|VTl- 

A summary of the fit results obtained for the pro- 
duction and rate separately for the neutrino, antineutrino 



and the combined neutrino and antineutrino samples is 
presented in Tab. IVIIII relative to the MC beam predic- 
tion. Fig.[T2]reports the correlation matrices respectively 
for the neutrino and the antineutrino muon angular dis- 
tributions. The full systematic covariance matrices for 
the neutrino and antineutrino angular distribution con- 
tain normalization uncertainties of 19% and 25%, respec- 
tively. For the initial cross section fits with pull terms, 
the common normalization uncertainty in the parameters 
is 6.5%, but as described above, the final fit for the if + 
weight determination only uses the parameter values and 
not the uncertainties. 

The average energy and angle for the com- 

bined neutrino and antineutrino samples are reported in 
Tab.llXl 

TABLE VIII. fit results for the rate and production rela- 
tive to the MC beam prediction for the neutrino, antineutrino 
and combined neutrino and antineutrino samples including 
the final x^/dof obtained from the K'^ production fit for NU- 
ANCE. Errors include statistical and systematic errors. The 
neutrino cross-section normalizations are held at the mini- 
mized values as listed in Table IVIII and are relative to the 
NUANCE predictions. 









Combined 




v-mode 


i^-mode 


u+i? mode 


K+ Prod. 


0.89±0.04±0.12 


0.54±0.09±0.32 


0.85±0.04±0.11 


K+ Rate 


0.94±0.05±0.11 


0.54±0.09±0.30 


0.88±0.04±0.10 


xVdof 


47.8/45 


18.5/27 


67.3/79 



TABLE IX. Measured -; — —, mean energy, and mean angle 
dpdil 

(with respect to proton beam direction) for the selected K'^ 
in neutrino, antineutrino, and the combined neutrino and an- 
tineutrino samples using NUANCE. Errors on the mean en- 
ergy and mean angle values correspond to the error on the 
mean for the relative distributions. 





Ek+ (GeV) 


Sk+ (degree) 


dpdO, 
(mb/(GeV/c x sr) 


v-mode 


3.84±0.03 


4.07±0.06 


5.77±0.83 


v-mode 


4.32±0.07 


2.04±0.09 


3.18±1.94 


V + v-madc 


3.93±0.03 


3.71±0.04 


5.34±0.76 



The fit results for the production double differ- 
ential cross-section in neutrino mode and antineutrino 
mode, though consistent, do not have to agree on the 
same central value. The sample of events selected 
in antineutrino mode has higher energy and lower angle 
with respect to the ones selected in neutrino mode as 
shown in Tab. IIVI and in Tab. IVIl Combining the neu- 
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Reconstructed Muon Angle Bin # 

(a) Neutrino Sample 




Reconstructed Muon Angle Bin # 

(b) Antineutrino Sample 

FIG. 12. Correlation matrix associated with the systematic 
and statistical uncertainties for the neutrino and antineutrino 
angular distributions using NUANCE. 



trino and antineutrino data gives a broader kinematical 
region for the rate and production measurements. 

The values for — — — for the neutrino, antineutrino, 
dpdi I 

and combined mode results are given in Tab. IIXI along 
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with the mean energy and angles for the corresponding 
K'^ samples. These values are obtained multiplying the 
measured K'^ production in Tab. IVIIII by the MC beam 
prediction in Eq. [TJ 

The reconstructed muon angular distribution for the 
SciBar 1, 2 and 3-Track sample rescaled using the fit 
results are shown in Fig. [T3] for neutrino mode and in 
Fig. HH for antineutrino mode. 




Reconstructed Muon Angle (degrees) 



(a) 1-Track Sample 




Reconstructed Muon Angle (degrees) 



(b) 2- Track Sample 




Reconstructed Muon Angle (degrees) 

(c) 3-Track Sample 




Reconstructed Muon Angle (degrees) 

(a) 1-Track Sample 




Reconstructed Muon Angle (degrees) 

(b) 2- Track Sample 




Reconstructed Muon Angle (degrees) 

(c) 3-Track Sample 

FIG. 14. Reconstructed muon angle for the SciBar 1-Track, 2- 
Track, and 3-Track samples after the fit in antineutrino mode 
running for NUANCE. The K'^ production weight and the 
cross-section values in Table IVIII have been applied to the 
NUANCE MC. The grey area represents the total systematic 
uncertainty in the MC. 



FIG. 13. Reconstructed muon angle for the SciBar 1-Track, 
2-Track, and 3-Track samples after the fit in neutrino mode 
running for NUANCE. The production weight and the 
cross-section values in Table IVIII have been applied to the 
NUANCE MC predictions. The grey area represents the total 
systematic uncertainty in the MC. 



VII. CONCLUSION 

In summary, we have made a measurement of the 
kaon production rate and cross section for 8 GeV pro- 



tons on a Be target in the Fermilab BNB using high- 
energy muon neutrino events observed in the SciBooNE 
detector. SciBooNE's full neutrino (antineutrino) data 
set, corresponding to 0.99 x lO^" (1.51 x 10^") POT, is 
used. The primary measurement uses the NUANCE in- 
teractions simulation and can be directly applied to the 
MiniBooNE oscillation analysis ^] which used that sim- 
ulation. A comparison of the results obtained with the 
NEUT neutrino interaction simulation gives similar re- 
sults (see Appendix [B|) . 

The analysis of the neutrino and antineutrino mode 
data are summarized in Tab. IVIIII Performing the 
analysis on the neutrino and antineutrino mode data 
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(Fa 

sets yields the following measurement for - — — for the 

dpdil 

p + Be — > + X at the mean energy of 3.9 GeV 
and angle 3.7 degrees: 



dpdVt 



(5.34 ± 0.76) mbl(GeV/c x sr) 



(5) 



The same analysis is also performed using the NEUT sim- 
ulation giving the results reported in Tab. IXVI and the 

resuh is: — ^ = {bA% ±0.1Q) mb/ {GeV / c x sr). As 
dpdil 

seen, the two simulations give very similar values (within 
2.5%) and demonstrate the independency of the result 
with respect to the underlying neutrino cross-sections 
and nuclear models that are different between NUANCE 
and NEUT. 

Applying these measurements will significantly reduce 
the systematic uncertainty associated with the measured 
Vf. background in MiniBooNE, which is a major source of 
uncertainty in the previously published oscillation ap- 
pearance result [Bj where the MiniBooNE collaboration 
used a conservative 40% error for the flux uncer- 
tainty. The rates (production x neutrino cross sec- 
tion) obtained using this analysis are: 

K+ Rate = 0.85 ±0.11, 

for both NUANCE and NEUT. In addition, the fact 
that the normalizations of the production and rate 
weights are less than one a away from the MC predictions 
(Eq. [T]) validates the procedure used to extrapolate the 
high energy kaon data to lower energy using the Feynman 
Scaling model [6] . The addition of this measurement in 
future fits to this model will reduce the systematic un- 
certainty associated with /v + production in low energy 
neutrino beams. 
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Appendix A: Cross-Section Fit Details 

In this section we will describe the details of the cross- 
section fit used in the minimization as described in 
Sec. Ell Eq. m The cross-sections described in Sec. IIIDI 
have been measured in the past few years with high pre- 
cision by the MiniBooNE coUaboration (2!. Isj-liH . 

Eq. lAll and Eq. IA2I are used to predict the number 
of neutrino events in bin i in case the fit is performed 
to measure the K'^ production or the K'^ rate. Slightly 
different equations, Eq. IA3I and Eq. IA4I are used in the 



antineutrino mode fit. i?]^™''' a-nd are the fit pa- 

rameters for the production and rate. 

We use different normalization factors depending on 
the neutrino parent and on the interaction type. The nor- 
malization factors depend on the neutrino parent {K^ ^ 
7r+, 7r~) and on the interaction types (QE, Itt, and 
Other). For production, the i?^^ and R],'" is ap- 
plied to both the and 7r+ while for rate, the 
R'^^ and R]^'^ is applied to only the 7r+ (the neutrino 
cross-sections are implicit in the normalization it- 
self). The reason we use the same normalization factors 
for neutrinos coming from tt"*" and is because we as- 
sume no dependence in the cross-section at high energy 
(> 2 GeV) from the neutrino parent. We apply different 
normalization factors between neutrino and antineutrino 
mode running for neutrinos from CCQE and CCItt in- 
teraction types to correctly take into account the Mini- 
BooNE measurement [4l| on the different neutrino flux 
prediction in antineutrino mode running. We use 
and R]j^ for neutrino mode running and i?^^^ and R^^"" 
for antineutrino mode running. We do not apply any 
rescaling to all the other neutrino interaction types. For 
neutrino mode running, we do not rescale any interac- 
tions coming from antineutrinos due to the very small 
contamination of antineutrinos. 



Prod. ^ 



f.Tpred\ TfRate 
JV \Rate = tC 



^+ x(7v:-'^-' + iyr--Kiv, 

7r+ 



if? 



i??^ X + X TV" 



(A2) 



For antineutrino mode running, we applied additional 
normalization factors to take into account neutrinos pro- 
duced by 7r~ as shown in Eq. IA3I and Eq. IA41 
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'QE 
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Ri^ X m: 



other 



P 



(A3) 



M^'-^'^lfl^.te = " X (A/; 



K 



R 



'QE 



X m;«^ 



Ry X A/, 



«^ + Af, 

+ 



"-4 



R^^ X A/p'^^ 



i?^'' X M, 



1", 



Ml 



other 



(A4) 



Additional pull terms are added to the fit when we 
fit the K'^ normalization factors and the neutrino cross- 
sections. We add a pull term for each cross-section weight 
(6 terms total) . These pull terms are then removed when 
performing the final minimization, and the cross sec- 
tion values are fixed at their best fit values. Each of the 



pull terms is of the form 



{Res - tJ.f 



where -Res is the 



cross-section best fit value, /i is the cross-section initial 
value, and a is the cross-section uncertainty. 

In the case of the antineutrino mode fit for the NU- 
ANCE MC we use a larger R^'^^^ and Rl'" with respect 
to the central fit values to get a better agreement be- 
tween data and MC for the SciBar 1-Track sample as 
shown in Fig. [9l For all other pull terms, ^ is set at 
1.0 with a 30% uncertainty. A summary of the cross- 
section best fit values is in Tab. En] for NUANCE MC 
and in Tab.|XlT]for NEUT MC. The variations in these 
normalization factors for the neutrino and antineutrino 
cross-sections are well within the cross-section systematic 
uncertainties and the cross-section best fit values agree 
well with the cross-section measurements performed by 
the MiniBooNE collaboration IsgMill. 




Reconstructed Muon Angle (degrees) 



(a) 1-Track Sample 




Reconstructed Muon Angle (degrees) 

(b) 2- Track Sample 




Reconstructed Muon Angle (degrees) 



(c) 3- Track Sample 

FIG. 15. Reconstructed muon angle for the SciBar 1-Track, 
2- Track, and 3-Track samples in neutrino mode running for 
data and NEUT MC. The background contributions from K~ 
and K1 are very small but included in the MC histogram. The 
grey area represents the total systematic uncertainty in the 
MC. 



Appendix B: K Measurements using NEUT 

In the SciBooNE experiment, neutrino interactions 
with carbon and hydrogen in the SciBar detector are also 
simulated by the NEUT [l3, [3 program library. 

The measurements presented in Sec. IIIII and Sec. IIVI 
are computed with assumptions on the underlying neu- 
trino cross sections and nuclear models. Due to that, 
we have repeated our analysis using NEUT [13, [13] as 
neutrino interactions Monte Carlo. The analysis is per- 
formed in the same way as for the NUANCE MC. The 
reconstructed muon angle distributions for the SciBar 1, 
2 and 3-Track samples are shown in Fig. \W\ for the neu- 
trino mode analysis and in Fig. for the antineutrino 
mode analysis. 



NEUT and NUANCE parameters are summarized in 
Tab. 1x1 The different parameter values result in different 
neutrino cross section predictions between the two and 
the differences are more in how the two Monte Carlos 
implement nuclear and cross section models. The total 
interaction numbers predicted by NEUT are summarized 
in Tab. |Xll 

The difference in the parameter choices fTab. lXj) . leads 
to NEUT predicting a larger QE and single pion rate 
than NUANCE. This difference in QE rate comes largely 
from the choice of k values, and the difference in single 
pion rate can be largely accounted for by the difference in 
M)^ assumptions. CC-QE interactions, which are imple- 
mented using the Smith and Moniz model [2l| . are the 
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TABLE X. Parameter used for Neutrino interaction simula- 
tion. 



(a) 1-Track Sample 




Reconstructed Muon Angle (degrees) 

(b) 2- Track Sample 




Reconstructed Muon Angle (degrees) 

(c) 3- Track Sample 

FIG. 16. Reconstructed muon angle for the SciBar 1-Track, 
2-Track, and 3- Track samples in antineutrino mode running 
for data and NEUT MC. The background contributions from 
K~ and K1 are small but included in the MC histogram. The 
grey area represents the total systematic uncertainty in the 
MC. 



dominant interaction in the SeiBooNE neutrino energy 
range. The nucleons are treated as quasi-free particles 
and the Fermi motion of nucleons along with the Pauli 
exclusion principle is taken into account. The Fermi sur- 
face momentum {pp) for carbon is set to 217(220) MeV/c 
and the nuclear potential (Eb) is set to 25(34) MeV/c for 
NEUT(NUANCE), as extracted from electron scatter- 
ing data [225. The default binding energy in NUANCE 
is somewhat higher than in NEUT because it addition- 
ally accounts for neutrino vs. electron scattering differ- 
ences [2^. With regards to the vector form factor, NEUT 
uses a dipole form with a vector mass of 0.84 GeV/c^, 
while NUANCE uses the BBA-2003 form factor [24|. A 
dipole form is used for the axial form factor with an 



Parameter 


NEUT 


NUANCE 


Pf 


217 MeV 


220 MeV 


Eb 


25 MeV 


34 MeV 


M^^-C 


1.21 GeV 


1.234 GeV 




1.0 GeV 


1.0 GeV 


K 


1.00 


1.022 


Mi- 


1.21 GeV 


1.10 GeV 


■^coherent 


1.03 GeV 


1.03 GeV 




(DIS) 


1.3 GeV 



TABLE XL The expected number of events in each neutrino 
interaction estimated by NEUT at the SciBooNE detector lo- 
cation with the neutrino beam exposure of 0.99 x lO'^" protons 
on target for neutrino mode and of 1.51 x lO'^'' for antineutrino 
mode. The 9.8 ton fiducial volume of the SciBar detector is 
assumed. Charged current and neutral current interactions 
are abbreviated as CC and NC, respectively. 



Interaction 



NEUT 

# Events # Events 

neutrino mode antineutrino mode 



CC QE 
CC single- TT 

CC Coh. TT 

CC DIS-t-Other 
NC 



50,330 
27,836 
1,692 
6,445 
35,153 



18,804 
7,350 
1,308 
1,760 
13,194 



adjustable axial mass, M^^ , for both NEUT and NU- 
ANCE. The values of M^^ = 1.21 GeV/c^ and k = 1.000 
(i.e. no additional Pauli blocking adjustment) are used 
in NEUT, and M^^ = 1.234 GeV/c^ and k = 1.022 
are used in NUANCE [2^. The same Fermi momentum 
distribution and nuclear potential are used in all other 
neutrino-nucleus interactions except for resonant tt pro- 
duction. 



The resonant production of single pion, kaon, and eta 
mesons, as described by the model of Rein and Sehgal 
(RS) [26|, is implemented with axial- vector form factors 
formahzed to be dipole with — 1.21 GeV/c^ for 

NEUT and M^" = 1.10 GeV/c^ for NUANCE. Reso- 
nance decays leading to multi-pion final states are also 
included in the NUANCE model as described in Sec. Ill Dl 
In NEUT, multi-pion production is simulated as deep in- 
elastic scattering (DIS) and the RS model is not used. 
Multi-hadron final states are simulated by two models: a 
custom-made program [i^ for events with W between 1.3 
and 2.0 GeV/c2 and PYTHIA/JETSET [H for events 
with W larger than 2 GeV/c^ The pion multiplic- 
ity is additionally restricted to be greater than one for 
1.3 < < 2 GeV/c^ to avoid double-counting sources of 
single pion production. In NUANCE, the DIS contribu- 
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tion slowly increases for W values starting at 1.7 GeV and 
becomes the only source of neutrino interactions above 
W >2 GeV. 

The total and inelastic pion-nucleon cross sections are 
treated differently in NEUT and NUANCE. In NEUT 
they are from the original Rein and Sehgal publication 
[H [IPl while in NUANCE, they are obtained from fits 
to PDG data 44,] and implemented as a function of pion 
energy with an additional rescale of the NC and CC co- 
herent pion production cross section. 

Intra nuclear interactions are simulated differently as 
well: more details can be found for both NEUT and NU- 
ANCE in [lllii. 

The results obtained for the cross-section best fit values 
using NEUT in the neutrino and antineutrino mode are 
summarized in Tab. IXIII 

TABLE XII. Summary of the fit results for the cross-section 
normalization factors as described in Appendix[X]with respect 
to NEUT predictions. The top four values are cross section 
normalization values for and i?^ coming from 7r+ and K'^ 
while the last two are normalization factors for coming 
from 7r~. 



TT+/K + 




Fit Value 




CCQE 


in 


f 


mode 


0.92±0.13 




CCQE 


in 


V 


mode 


1.05±0.25 




CClvr 


in 


V 


mode 


1.12±0.25 




CClTT 


in 


V 


mode 


0.94±0.27 










CCQE 


in 


V 


mode 


1.01±0.30 




CClTT 


in 


v 


mode 


1.01±0.30 



The mean energy and mean angle (with respect to pro- 
ton beam direction) for the selected and the mean 
energy for z^^ from the selected in each of the three 
samples, for both neutrino mode running and antineu- 
trino mode running, are summarized in Tab. IXIIII and 
Tab. IXIVI A summary of the fit results obtained for the 
K'^ production and rate separately for neutrino, antineu- 
trino and the combined neutrino and antineutrino sam- 
ples is presented in Tab. IXVI relative to the MC beam 
prediction. Fig. [17] reports the correlation matrices re- 
spectively for the neutrino and the antineutrino mnon 
angular distributions. For NEUT, the full systematic co- 
variance matrices for the neutrino and antineutrino an- 
gular distribution contain normalization uncertainties of 
19% and 25%, respectively. For the initial cross section 
fits with pull terms, the common normalization uncer- 
tainty in the parameters is 6.5%. 

The average K'^ energy and angle for the com- 
bined neutrino and antineutrino samples are reported in 
Tab. lXVn 

The values for — — — for the neutrino, antineutrino, 
apai I 

and combined mode results are given in Tab. IXVTl along 
with the mean energy and angles for the corresponding 



TABLE XIII. NEUT MC predicted mean energy and mean 
angle (with respect to the proton beam direction) for the se- 
lected -K"*" samples and predicted MC energy for the selected 
from in the three SciBar samples in neutrino mode 
running. Errors correspond to the error on the mean values. 
The RMS of the relative distributions is reported in square 
parenthesis. 





£;K+[RMS](GeV) 


[RMS](deg) 


i;.,JRMS](GeV) 


1-Trk 


3.6±0.1[1.2] 


4.2±0.1[2.0] 


3.1±0.1[1.0] 


2-Trk 


3.9±0.1[1.2] 


4.0±0.1[2.0] 


3.3±0.1[1.0] 


3-Trk 


4.2±0.1[1.1] 


3.6±0.1[1.9] 


3.6±0.1[0.9] 



TABLE XIV. NEUT MC predicted mean energy and mean 
angle (with respect to the proton beam direction) for the se- 
lected K'^ samples and predicted MC energy for the selected 
from for the three SciBar samples in antineutrino mode 
running. Errors correspond to the error on the mean values. 
The RMS of the relative distributions is reported in square 
parenthesis. 





£K+[RMS](GeV) 


ej^+[RMS](deg) 


£.,JRMS](GeV) 


1-Trk 


4.1±0.1[1.3] 


2.2±0.2[1.6] 


3.1±0.1[1.1] 


2-Trk 


4.3±0.1[1.2] 


1.7±0.1[1.2] 


3.4±0.1[1.0] 


3-Trk 


4.7±0.2[1.1] 


1.6±0.2[1.0] 


3.8±0.2[1.0] 



TABLE XV. fit results for the rate and production rel- 
ative to the MC beam prediction for the neutrino, antineu- 
trino and combined neutrino and antineutrino samples includ- 
ing the final x^/dof obtained from the production fit for 
NEUT. Errors are statistical and systematic errors. The neu- 
trino cross-section normalizations are held at the minimized 
values as listed in Table IXIII and are relative to the NEUT 
predictions. 



Combined 

i^-mode u-maAR v^u mode 

A'+ Prod. 0.90±0.05±0.13 0.77±0.12±0.31 0.87±0.05±0.11 
A'+ Rate 0.92±0.05±0.11 0.76±0.12±0.27 0.88±0.05±0.10 
xVdof 40.6/45 17.7/26 58.9/77 



samples. These values are obtained multiplying the 
measured production in Tab. IXVI by the MC beam 
prediction in Eq. [T] 

The reconstructed muon angle distributions for the 
SciBar 1, 2 and 3- Track samples rescaled using the fit 
results are shown in Fig. [18] for the neutrino mode anal- 
ysis and in Fig. [19] for the antineutrino mode analysis. 
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TABLE XVI. Measured 



(fa 



mean energy, and mean angle 



dpdQ 

(with respect to proton beam direction) for the selected K'^ 
in neutrino, antineutrino, and the combined neutrino and an- 
tineutrino samples using NEUT. Errors on the mean energy 
and mean angle values correspond to the error on the mean 
values of the relative distributions. 





Ek+ (GeV) 


Qk+ (degree) 


dpdQ. 
(mb/(GeV/c x sr) 


i^-mode 


3.85±0.03 


4.02±0.06 


5.85±0.90 


u-mode 


4.29±0.07 


1.93±0.09 


4.59±1.97 


V + v-maA& 


3.92±0.03 


3.66±0.04 


5.49±0.79 




5 10 15 20 25 30 35 
Reconstructed Muon Angle Bin # 



(b) Antineutrino Sample 



FIG. 17. Correlation matrix associated with the systematic 
and statistical uncertainties for the neutrino and antineutrino 
angular distributions for NEUT. 
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FIG. 18. Reconstructed muon angle for the SciBar 1-Track, 
2-Track, and 3- Track samples in neutrino mode running for 
data and the NEUT MC. The K+ production weight and the 
cross-section central values in Table IXIII have been applied 
to the NEUT predictions. The grey area represents the total 
systematic uncertainty in the MC. 



FIG. 19. Reconstructed muon angle for the SciBar 1-Track, 
2-Track, and 3- Track samples in antineutrino mode running 
for data and NEUT MC. The K+ production weight and the 
cross-section central values in Table IXIII have been applied 
to the NEUT predictions. The grey area represents the total 
systematic uncertainty in the MC. 
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